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Abstract 

We present results from a numerical study of AA = 1 supersymmetric Yang-Mills theory using do- 
main wall fermions. In this particular lattice formulation of the theory, supersymmetry is expected 
to emerge accidentally in the continuum and chiral limits without any fine-tuning of operators. 
Dynamical simulations were performed for the gauge group SU{2) on 8^ x 8 and 16'^ x 32 lattice 
space-time volumes and at three different values of the coupling: f3 = 2.3, 2.353 and 2.4. Results 
from this study include measurements of the static potential, residual mass, and a chirally extrap- 
olated value for the gluino condensate at /? = 2.3. In addition to these, we study the low lying 
eigenvalues and eigenvectors of the five dimensional Hermitian domain-wall fermion Dirac oper- 
ator and present evidence that, for the choice of parameters under investigation, features of the 
spectrum appear qualitatively consistent with strong coupling and the presence of a large residual 
mass. From the five dimensional eigenvalues we explore the possibility of using the Banks-Casher 
relation to determine an independent value for the gluino condensate in the chiral limit. 
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I. INTRODUCTION 



In recent years, substantial effort has been devoted toward formulating supersymmetric 
(SUSY) gauge theories on the lattice.^ These efforts have been partially motivated by the 
fascinating theoretical and technical difficulties associated with the problem, as well as the 
obvious potential role of SUSY in beyond the standard model physics. Of crucial importance 
pertaining to the latter point is an understanding of a variety of nonperturbative phenomena, 
including-but not limited to-dynamical SUSY breaking. An understanding of nonperturba- 
tive aspects such as these may in principle be achieved with numerical simulations, provided 
an appropriate lattice discretization of the theory may be found. 

Since naive lattice discretizations typically break SUSY explicitly, simulations of such 
theories will generally require a substantial degree of fine-tuning in order to cancel any 
undesirable SUSY breaking operators which may arise through radiative corrections. Such 
a task is exceedingly difficult, even when the number of parameters which require fine- 
tuning are relatively small. However, it has been realized for some time that one of the 
simplest SUSY theories, M = 1 supersymmetric Yang-Mills (SYM), may be simulated using 
conventional lattice discretizations and yet require only a minimal degree of fine-tuning. 

The field content of A/" = 1 SYM consists of a vector field and a single Majorana fermion 
which transforms as an adjoint under the gauge group, which for our purposes will be 
taken to be SU{N^. The theory possesses an anomalous U{1)a axial symmetry, however, 
a discrete Z2Nc subgroup of U{1)a survives at the quantum level. It is believed that N = 1 
SYM shares a var.ety of features i,r wth Quantum Chromodynamies (QCD). such 
as confinement and chiral symmetry breaking [3]. But in contrast to QCD, A/ = 1 SYM 
is believed to possess discrete chiral symmetry breaking from N2Nc —>■ Z2 which results in 
the possible formation of domain walls but no Goldstone bosons. Although this theory does 
not possess dynamical SUSY breaking, as implied by a nonvanishing Witten index ^, it 
is still believed to exhibit a variety of other features which may be of interest to explore 
nonperturbatively. Of perhaps greatest interest is whether or not a gluino condensate forms, 
the spectrum of the theory, and the relationship between the domain wall tension and the 
gluino condensate. 



^ For recent reviews, see e.g., P, 0]. 
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The low energy spectrum of = 1 SYM is believed to consist of massive supermultiplets 
which may involve glue-glue, glue-gluino as well as gluino-gluino bound states. Although 
the states within a given supermultiplet are degenerate, at finite gluino mass one expects 
mass splittings which may be calculated using a variety of effective theories 5|, ISj. Lattice 
simulations of A/" = 1 SYM may be used to directly test the predictions of these effective 
theories, as well as to determine the effects of soft SUSY breaking on observables other than 
the spectrum. 

N = 1 SYM theory provides an ideal starting point for studying aspects of SUSY numer- 
ically while still allowing one to utilize familiar and well understood fermion discretizations. 
Because the field content of this theory consists only of a vector field and Majorana fermion 
(and in particular no scalar fields), with conventional lattice discretizations of this theory, 
the only relevant SUSY violating operator which may arise radiatively is a gluino mass term. 
In the chiral and continuum limits, SUSY is therefore expected to emerge accidentally at 
infinite volume. For fermion discretizations that explicitly break chiral symmetry, taking 
the chiral limit requires tuning the input mass parameter to some critical value, such that 
the bare mass vanishes. 

n n 

With the use of a Ginsparg- Wilson ['3] fermion discretization such as overlap fermions [8| , 
fine-tuning of the gluino mass term may be entirely avoided. This study utilizes the closely 



which realizes chiral symmetry 



related domain wall fermion (DWF) formalism 

ay 

in the limit of infinite fifth dimension extent and vanishing input gluino mass. In a typical 
DWF simulation where the fifth dimension is finite, however, the action receives residual 
chiral symmetry violating radiative corrections, including an additive mass renormalization. 
The important advantage that DWF fermions have over other discretizations (which lack 
chiral symmetry) is that these renormalizations may be removed by taking a controlled and 
unambiguously determined limit. 



In the past, a variety of numerical studies have emp 



simulate U = 1 SYM 
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This discretization requires fine-tuning and formally possesses a sign problem, since the 
Pfafiian obtained from integrating out the fermion degrees of freedom is generally not positive 
definite.^ Studies, however, have found that the Pfaffian phase becomes innocuous in the 



^ While one may prove that the fermion Pfaffian for this theory is positive definite in the continuum, 
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parameter regime of physical interest and may be effectively handled by phase reweighting 



techniques. It was observed in 



301 1 that DWFs have an additional advantage over Wilson 



fermions in that the fermion Pfaffian is positive definite. With the use of DWFs, one may 
therefore avoid the task of phase reweighting altogether. 

The first and until recently the only study to use DWFs to investigate M = 1 SYM 
focused on the chiral limit of the gluino condensate [311] at a single lattice spacing. This 
study was performed by using an inexact, hybrid molecular dynamics R (HMDR) algorithm 



32l | which, for the theory under consideration, is prone to finite integration step size errors. 



Our study offers an im pro vement over 3l|], by utilizing the rational hybrid Monte Carlo 
(RHMC) algorithm 33|, |34, l35| , which is an exact algorithm. Furthermore, due to algorith- 



mic improvements and faster computers, we are able to explore the theory on larger lattices 
and smaller couplings as well. Recently an independent study of A/" = 1 SYM using DWFs 
was reported by Giedt, et al. in 



36 



37|. 



The primary purpose of our study is to establish a set of sensible simulation parameters, 
perform basic measurements and provide the necessary ground work for more detailed studies 
of the theory in the future. As such, we expand the work of 3l| in several important respects: 
we 1) establish the lattice scale by measuring the static potential and provide evidence for 
confinement which is consistent with expectations, 2) determine the size of the residual mass 
in order to ascertain the proximity to the SUSY point, 3) extrapolate the chiral condensate 
to the chiral limit using a recent, theoretically motivated fit formula for its Ls dependence 



and which differs from the formula used in 



31| . 4) study the structure of the low lying 



eigenvectors and eigenvalues of the Hermitian DWF Dirac operator and attempt to extract 
an independent value for the chiral condensate using the Banks-Casher relation 



. With 



exception to the third and final points, questions such as these could not be easily addressed 



m 



3l| due to the limited space-time volumes employed. 



Portions of this work have been reported at Lattice 2008 [39|, including results obtained 
for points 1-3. In this paper we go into greater detail on the analysis of these results, and 
expand this analysis by including the eigenvalue and eigenvector studies outlined in point 4. 
We organize this paper as follows. In Sec. [Ill we describe the details of our numerical simu- 
lations, including a brief review of the domain wall fermion formalism and conventions, the 



positivity is not necessarily guaranteed on the lattice. 
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RHMC algorithm, ensembles and measurement methods. Sec. IIIII is devoted to issues con- 
cerning thermalization and autocorrelations, while Sec. IIVI pertains to measurement details, 
data analysis and results. In Sec. |Vl we summarize our findings. Finally, in Appendix |X] 
we review some basic properties of Marjorana fermions, which may be less familiar to those 
readers who specialize in lattice QCD, and specifically derive expressions for some of the 
gluino-dependent observables studied in this paper. 

II. SIMULATION DETAILS 

A. Lattice action and conventions 

Dynamical numerical simulations of A/" = 1 SYM were performed using a Wilson gauge 
action and domain wall fermions. The full form of the partition function in Euclidean 
space-time is given by: 



where Sq represents the gauge action and Sp represents the action for domain wall fermions. 
Spv represents the action for the Pauli-Villars (PV) fields which are required in order to 
cancel off the UV contributions associated with the bulk fermions in the fifth dimension. 
Simulations are performed on a x lattice, where V = x T represents the lattice 
space-time volume of the four physical dimensions with spatial extent L and temporal extent 
T, and Lg represents the size of the unphysical fifth dimension. 
The Wilson gauge action is given by: 



where Px,fiu = Tr {Ux^^iUx+^^yUlj^^^^U}.^^) represents a 1 x 1 plaquette and Ux,^ are link vari- 
ables which belong to fundamental representation of SU{Nc) and are associated with the 
four dimensional space-time coordinate (x) and orientation (/i). As is usual with the domain 
wall fermion discretization, the gauge field is taken to be independent of the fifth dimension 
coordinate (s). The bare coupling is given by /? = i/g"^. 

The five dimensional DWF Dirac operator D is defined by: 




(1) 




(2) 



D 



(3) 
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where 



—mfP_ , if s' = and s = Lg — 1 

-nifP^ , if s = and s' = - 1 (5) 

P-6s+i,s' + P+Ss-i,s' - 26s,s' , otherwise , 



and the spin projectors P±^ and are given by: 

^± = ^(1±75), P^ = lil±l,)- (6) 

The gamma matrices 7^ in Euchdean space are taken to be Hermitian and satisfy the proper- 
ties: Tr (7^71/) = 4:6^1, and (7^^)^ = 1 (no sum on /i); the four dimensional chirahty operator 
is given by 75 = 70717273- The real link variables V.,^ belong to the adjoint representation 
of the gauge group. These may be expressed in terms of fundamental representation link 
variables t/.,^ which appear in the gauge action via the identity: 

yab ^ 2Tr {T^U^T'^U) , (7) 

where are generators of the fundamental representation of the gauge group which satisfy 
Tr (T"T^) = 5'^V2, with a = 1, . . . , A^^ - 1. Finally, M5 and rrif appearing in Eq. [3] represent 
the domain wall height and input gluino mass, respectively. 
The DWF Dirac operator satisfies the following properties: 

rlDT, = , 

C^DC = D* , (8) 

where Fs = -R575, -R5 is the reflection operator about the mid-plane in the fifth direction, 
C = R^C is the DWF charge conjugation operator and C is the standard four dimensional 
charge conjugation operator which satisfies the usual relations: 



C^l.C = -7j , 
C^7.C = 75^ , 



P± — Pr/l in the notation of, for example, reference [4o| . 



C^C = 1 , 

c^t:,^c = -n^r . (9) 

Here, T^^^ are the generators of the adjoint representation of the gauge group which satisfy 
TTT;^ajTl^aj = N^d^K The DWF and PV actions are given by: 

Sp[^, U] = ^^D(M5, mf)m , Spv[<l>, U] = hD{M,, 1)$ , (10) 

where we have imposed the DWF and PV Majorana conditions: \E' = "^^C and $ = <I>^C 
Note that the PV action used in these simulations is not the one introduced in kOl but 
rather a variant of that action [41]. 

All gluino dependent observables in this study are expressed in terms of the four di- 
mensional gluino interpolating fields q{x), which are defined at the boundaries of the fifth 
dimension. In terms of the five dimensional fields \E'(a;, s), the gluino interpolating fields are 
given by: 

q{x) = P_^(x,0) + P+^(x,L, - 1) , 

q{x) = ^(x, L, - 1)P_ + ^(x, 0)P+ . (11) 

For convenience, we also define here the four-dimensional gluino fields Qmi^) which may be 
associated with the mid-plane in the fifth dimension: 

qm{x) = P_^(x,L,/2) + P+^(a;,L,/2-l) , 

q^x) = ^(x, Ls/2 - 1)P_ + ^(x, LJ2)P+ . (12) 

Note that each of the interpolating fields in Eq. [11] and Eq. [12] satisfy the appropriate four 
dimensional Majorana conditions: q = q^C and g^, = q^C. 



Following 

by: 



401 1 we may define four- dimensional vector and axial currents V^(x) and A^{x) 



Ls-l 



s=0 

= XI sgn U - ^^^y^ ) jM(x, s) , (13) 

where j^(x, s) are the first four components of the five dimensional conserved current, given 
by: 

j^(x, s) = ^(x + fi, s)P+l^,^^^(x, s) - ^(x, s)p-V^,^^{x + fi,s) , (14) 



and the last component is given by: 

{^(x, s)P_^(x, s + 1) - ^(x, s + l)P+^(x, s) , < s < - 1 , 

_ _ (15) 

^(x, Ls - l)P_^(x, 0) - ^(x, 0)P+^(x, Ls-l) , s = Ls-l . 

The divergence of the four dimensional vector and axial currents satisfy 

A^V^(x) = , 

A^^^(x) = 2m/J5(a;) +2J5q(x) , (16) 
where A^/(x) = /(x + /i) — /(x) is the forward finite difference operator and 

M^) = 3b{.x,Ls - 1) = g(x)75g(x) , 

J5q{x) = j5(x,L,/2 - 1) = qmix)'y5qm{x) , (17) 

are the pseudo-scalar densities defined on the walls and midpoint of the fifth dimension, 
respectively. The axial Takahashi-Ward identity (ATWI) is given by: 

A^{A,0) = 2mf{J,0) + 2(J5,0) + i{50) , (18) 

where it may be verified that the term involving the midpoint pseudo-scalar density J5q{x) 



gives rise to the desired axial anomaly 42|. 



B. Simulation method 

In order to simulate the partition function given by Eq. [H it is necessary to first integrate 
out the DWF and PV degrees of freedom. After performing this integration, we are left with 
an effective action which depends only on the gauge fields and is given by: 

Pf[CD{M„l)] ' ^'^^ 
where, using the relation detC = 1, the Pfaffian may be expressed as 

Pf (CD) = VdetD . (20) 

Numerical simulations of = 1 SYM requires a positive definite fermion Pfaffian in order 
to unambiguously define the Pfaffian as the square root of a determinant. It may be shown 
that for domain wall fermions, this is indeed the case and demonstrates an advantage of 

8 



using this formalism over other discretizations, such as Wilson fermions 30|. By exploiting 
Fs-Hermiticity described by Eq. [8], we may rewrite the effective action as 



1/4 



(21) 



where T){M^^mf) = D{M^^mf)'^ D{M^^mj). For the purpose of numerical simulation, we 
introduce a single complex pseudo-fermion field x order to reproduce the effects of the 
ratio of fermion determinants given by Eq. EH"^ The resulting partition function is given by: 

[dU] [rfxl [dx] e-^«[^]-^^^[^'''^'^] , (22) 

where the pseudo-fermion action Spp is given by: 

Spf[x\xM = X^V'/\M,,l)V-'/\M,,mf)V'/\M,,l)x • (23) 

Dynamical numerical simulations of Eq. [22] with Nc = 2 were performed using a modified 
version of the Columbia Physics System (CPS)-a software system which is developed and 
maintained by the RBC collaboration for the purpose of studying lattice Quantum Chro- 
modynamics. Modifications to the software were specifically made in order to accommodate 
the adjoint character of the fermions being simulated as well as to reduce the gauge group 
under consideration from SU{3) down to SU{2). We perform numerical simulations of the 
partition function given by Eq. [22] via the exact rational hybrid Monte Carlo (RHMC) algo- 



rithm 



33 



34 



35| . Details of the parameters used in the approximation of V^{M^,mf) and 
©^(Ms, 1), for the appropriate rational powers p and q (i.e., p = 1/8 and g = 1/8 for the 
pseudo-fermion refreshment step, and p = 1/4 and g = 1/8 for the evolution) used in this 
simulation, may be found in Table [T] and Table [Til Parameters attributed to the PV fields in 
this study are labeled with the superscript PV. The parameters Xmin (A^^) and Xmax (AmD 
specify the lower and upper bounds on the eigenvalue range of ^{M^^mf), over which we 
require the rational approximation to be valid. The rational approximation is used both in 
the molecular dynamics (MD) evolution of the gauge field as well as in the accept/reject 
Monte Carlo (MC) step, which is performed at the end of the trajectory in order to remove 
any errors associated with the finite step size 5t in the MD evolution. The parameters 



^ The simultaneous treatment of normal and PV fermion contributions with a single pseudo-fermion field 
was an idea of M. Clark, and was first introduced in |43|. 
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f^MD (jI'^d) umc iP'^o) represent the degree and therefore accuracy of the rational 
approximation in each of these steps. In our simulations we require a greater accuracy in 
the accept/reject step compared to the evolution and therefore take nMc > ^md- Similarly, 
for the MD evolution and MC accept/reject step, the conjugate gradient stopping conditions 
used in the operator inversions were 1 x 10~^ and 1 x 10^^°, respectively. 

Finally, the MD evolution was performed using a two level Omelyan integrator Q, 
with the Omelyan integrator parameter set to A = 0.215. The length of each trajectory in 
the evolution is given by 5 x 5r MD time units. At the end of each trajectory, we project 
the gauge field back onto the gauge group. This step breaks reversibility, however the errors 
induced by the projection are negligible. Parameters and characteristics of the MD evolution 
and accept/reject MC step are provided in Table UTTl and Table HVl In each ensemble, the 
time step 6t was chosen such that the acceptance rate was approximately 70 — 80%, however 
for some ensembles the acceptance rates were as low as 60% or as high as 90%. One may 
verify from Table UTTl and Table HVl that the equality {e~^^) = 1 holds within errors for all 
ensembles, which indicates that the algorithm is working correctly with our modifications 



to the code 



46|. 



C. Ensembles and simulation parameters 

Numerical simulations where performed on two different lattice space-time volumes: V = 
8^ X 8 and 16'^ x 32. For the gauge fields we impose periodic conditions (BCs) in all space- 
time directions, whereas for the gluino we impose periodic BCs in the spatial directions and 
anti-periodic BCs in the temporal direction. The 8^ x 8 ensembles were generated using the 
parameters: (3 = 2.3, Lg = 12, 16, 20 and 24, rrif = 0.02 and M = 1.9, in order to check our 
code against the results of ^] . Simulation parameters for these configurations are listed in 
Table [T] and Table IIIII Estimates of the gluino condensate and average plaquette are given 



in Table |Vl and are consistent with 31i at the l-2cr level. 

Our 16^ X 32 ensembles were generated at /? = 2.3, Lg = 16, 20, 24 and 28, and input 
gluino masses: ruf = 0.01, 0.02 and 0.04. As with the 8'^ x 8 ensembles, the domain wall 
height is set to M = 1.9. For j3 = 2.3 and ruf = 0.02, we extended the range of Lg values to 
include Lg = 32, 40 and 48 in order to investigate the dependence of the residual mass and 
gluino condensate on the size of the fifth dimension. Several simulations where performed 
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at the weaker couplings j3 = 2.353 and f3 = 2.4 as well in order to investigate the coupling 
dependence of the residual mass and eigenvalues of the five dimensional Hermitian DWF 
Dirac operator. Table [III and Table [IV] provide a complete list of the ensembles associated 
simulation parameters used. Gluino condensate and average plaquette results for these 
ensembles are provided in Table IVI[ 

D. Measurement methods 

1. Residual mass 

At long distances, finite lattice spacing effects may be characterized by a continuum 
Symanzik effective Lagrangian: 

^Symanzik = ^SYM + CL ^ ^-1 + O.C1 + . . . , (24) 

where the leading contribution Csym represents the target continuum Af = 1 SYM theory. 
The contributions with n = — 1, 1, . . . characterize finite lattice spacing effects at order 
C(a"), up to possible logarithms. Specifically, the lowest order contributions to Csymanzik 
may be expressed as: 

a"^£_i = ^(m/ + mres)J^ , aCi = ^c.^Js/ , (25) 

where 

^5 = , ^5/ = r''r^,.F'^,r , (26) 

and ip and F^^^, are the four dimensional continuum fermion field and color field strength 
tensor. The invariant tensor /"^'^ = — 2i Tr T'*[T'', T^] represents the structure constants of 
the gauge group. The residual mass rrires in Eq. [25] characterizes the leading chiral symmetry 
breaking effects due to the finite extent of Lg, and is defined in such a way that the bare 
gluino mass is given by the simple sum: 

rrig = rrif + rrires ■ (27) 

The Sheikholeslami-Wohlert term which is proportional to c^^ will in turn depend on both 
mf and L^, and is expected to vanish in the — and Lg ^ oo limits. 
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In the continuum effective theory, the ATWI will read: 

A^{A^{x)0{x')) = 2{mf + mres){Ux)0{x')) + {ptop{x)0{x')) + i{50{x')) + . . . , (28) 
where 

TV 

Aop = ^TrFF, (29) 

and the contribution proportional to Csw has been omitted as it is higher order in the lattice 
spacing and is suppressed at large L^. Comparing the continuum ATWI with that of the 
lattice expression, we expect at long distances and sufficiently close to the continuum limit 
the identity: 

J->q{x) niresMx) + Ptop{x) . (30) 

It should be emphasized that there are no anomalous contributions to the gluino mass (as 
there may be for the quark mass in one ffavor QCD) due to the underlying Z2Nc symmetry 
which is present in the target A/" = 1 SYM theory. 

Next we describe details on how the residual mass may be extracted from the the low 
energy identity Eq. [301 obtained above. The cleanest method for extracting the residual 
mass in the context of QCD is to study the long time behavior of the flavor non-singlet wall- 
midpoint pseudo-scalar density correlator divided by the wall-wall pseudo-scalar density 
correlator. In the case of A/" = 1 SYM, however, the method is complicated by the fact 
that there are no flavor non-singlet pseudo-scalars, and that the flavor singlet pseudo-scalar 
correlator is contaminated by anomalous contributions. 

In order to proceed, we consider separately the connected and disconnected contributions 
to the wall-midpoint and wall-wall pseudo-scalar density correlators, which may be defined 
in a standard fashion. For the wall-midpoint correlator we have: 

{J,,{x)P{x')) = 2 ( J5.(x)P(a:'))„,,, + {Mx)P{x')),„^,^, , (31) 

where, if the pseudo-scalar density P{x) is taken to be J5(x), for example, the connected 
and disconnected parts of the wall-midpoint correlator are given by: 

/c ^ ^ 

( (a;) Js (a;')) connected = ( 9m(a^)75 gm(a^)g(a^0^ 759(a^0 

(^5<?(x)J5(x'))di™ected = ( srnixh5qUx)q{x')-f5q{x'l ) , (32) 
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respectively. In each of these equations the braces indicate which gluino fields are contracted. 
The factor of two appearing in the connected contribution to Eq. [3T] above accounts for the 
fact that we are working with Majorana rather than Dirac fermions (i.e., the quark fields 
q and q may be contracted with themselves; for more details, see Appendix Similar 
expressions may be defined for the wall-wall pseudo-scalar correlator. 
At low energies and long distances we expect the relation: 

{J^g{x)P{x')) ^ mres{Mx)P{x')) + {ptop{x)P{x')) , (33) 

which we assume for the moment may be decomposed in terms of connected and disconnected 
parts as: 

{J^qi.x)P{x )) connected ~ '"^'"es (-^5 (a^)-P(a; ) )conjierfed ' 

disconnected 

+ {pto,{x)P{x')) . (34) 

If this decomposition is indeed valid, then the residual mass may be extracted from the ratio 
of long time correlation functions: 

mUmf) = lim R{t) , i?(t) = (Ex '^s.lx, t)P(0))„,^, ^^^^ 



(Ex^5(x,t)P(0))_,^, 



Note that there is an 0{a) ambiguity in this particular expression for the residual mass 
due to the presence oi nif and dependent higher order contributions to Eq. [28] which 
have been omitted in Eq. [301 As such, we adopt the notation m^gg(m/) rather than rrires to 
represent the quantity which is extracted from the ratio R{t). 

In the case of two- or two plus one-fiavor QCD, the correctness of the decomposition given 
in Eq. [3ll may be proved trivially because there is an underlying fiavor symmetry which 
allows one to relate the connected contributions to these diagrams to their non-anomalous 
fiavor non-singlet counterparts. In A/" = 1 SYM, however, there is no such underlying 
fiavor symmetry and therefore Eq. [35] remains based upon an assumption. One may verify 
from numerical simulations that the ratio R{t) indeed tends to a constant value, which 
suggests that the time dependence of the connected part of the wall-midpoint and wall-wall 
correlators are the same. However this observation does not preclude the possibility that 
the disconnected contribution to the wall-midpoint correlator on the left-hand side of Eq. [33] 
contributes to the connected part of the wall-wall correlators on the right-hand side at low 
energies, thus giving an incorrect estimate of the residual mass. Ultimately, establishing 
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the validity of Eq. [35] will require additional theoretical analysis or corroboration from an 
independent calculation of the residual mass. Close to the continuum limit one may, for 
instance, be able to compute the residual mass from the valence mass (m^) dependence of 
the Hermitian DWF Dirac operator eigenvalues and verify that the value extracted from 
R{t) measurements, modulo 0{a) corrections, is consistent. Details of this approach are 
explained in greater detail in Sec. IIID4I 



2. Static potential 

The potential associated with two static, fundamental representation sources separated 
by a distance |x| may be obtained from the Wilson loop: {W{'x,t)). Wilson loops were 
measured using Coulomb gauge fixed gauge field configurations with link fields belonging to 
the fundamental representation of the gauge group. In order to reduce the statistical errors 
associated with this observable, Wilson loops were measured with their time axis oriented 
along each of the four directions of the lattice and the results were then averaged. The 
space-time dependence of the Wilson loop is expected to be of the form: 

(Vr(x, t)) = C(x)e-^('')* + excited states , (36) 

and from this the static potential V^(x) may be extracted from: 

at asymptotically large times. The signal to noise ratio associated with Eq. [37] generally 
deteriorates in this limit, therefore moderate values of time are instead used to extract the 
potential. 

In this work, we use Eq. [37] to locate the plateau region of Ve//(x, t) for each value of x. 
We then fit the Wilson loops to the functional form given by Eq. [3S] in order to extract the 
potential. Finally, the extracted values of V^(x) may be fit to the Cornell potential which is 
given by: 

(y 

l^(x) = K)- — + a|x| , (38) 

allowing us to determine the constant term (Vq), Coulomb term (a), string tension (a) and 
Sommer scale (rg) defined by : 



d\x\ 



= 1.65 , (39) 

|x|=ro 
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at fixed values of tlie lattice spacing and gluino mass. 



3. Gluino condensate 
The gluino condensate 

X 

was measured using a stochastic estimator with a single random Gaussian space-time volume 



source. Note that we normalize the gluino condensate following the conventions of [48|] (i.e., 
we divide by 4 spin x 3 adjoint color components), such that (qq) ~ mj^ for large mass. For 
these measurements, the conjugate gradient stopping condition used for the inversion of the 
Dirac operator was set to 1 x 10^^^. Measurements of the condensate were made for several 
different values of m/ and Lg, then a chiral limit extrapolation of the data was performed. 
Details of this analysis are provided in Sec. IIV CI . 

A second, independent measurement of the chiral condensate may be obtained using the 



Banks-Casher relation 38|], which in the continuum reads: 



lim lim -{^^) = ^p(O) , (41) 

where is a. continuum Majorana gluino field, and p(A) is the density of eigenvalues per 
unit volume of the four-dimensional continuum Dirac operator^, given by: 

p(A) = lim hm i { 5^5(A - A.)\ . (42) 

In Sec. HID 41 we describe how one may obtain a numerical estimate of p(A) by studying the 
low-lying spectrum of the five dimensional Hermitian DWF Dirac operator. In this study we 
are unable to measure p(A) per se, but rather a closely related quantity p'(A;mg), which is 
given by Eq. |42]prior to taking the V ^ oo and nig limits. Note that p'(A; m^) implicitly 
depends on the gluino mass through the distribution of gauge fields which has been sampled 
at finite ruf and Lg. The gluino condensate at a finite mass may be extracted from p'(0; m^) 
in the infinite volume limit, in a way very much analogous to Eq. |4T], by taking the valence 
mass goes to zero limit. Unlike the condensate obtained from Eq. HO, however, the condensate 
extracted by this method is free of UV divergent contributions and, as such, should only 
depend on mj and Lg in the particular combination: nig = mj + nires, provided higher 
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order terms in the Symanzik action may be neglected. Hence, by measuring this quantity 
as opposed to Eq. HQl we ehminate having to perform two independent extrapolations (i.e., 
nif —' and Lg — > oo) in favor of a single extrapolation of nig 0. 



4. Eigenvalues of the Hermitian DWF Dirac operator 

Here we describe some properties of the Hermitian DWF Dirac operator and its eigen- 
values. Before proceeding, however, we first review some of the expected features of the 
continuum, four dimensional Dirac operator D'^ =11) + nig and the four dimensional Hermi- 
tian Dirac operator = Using the anti-Hermiticity oi p = 7/^-0/^, one may show 
that analogous to Eq. [HI D'^ satisfies the following relations: 

C^D^C = {D^ . (43) 

From these properties it is easy to show that the eigenvalues of come in complex conjugate 
pairs given by ±zA + m^, where ±zA are eigenvalues of ij) and A G M. Furthermore, the 
eigenvalues of are two- fold degenerate. The eigenvalues of D\j, which we shall denote 
\hi are real and given by ±a/A^ + ni^g- Similarly, the eigenvalues of also have a two-fold 
degeneracy which follows from the property 

{l,C)^ D%{^,C) = {Dj,r . (44) 

Specifically, if D^jipx^ = \h'>P\hi then it follows that D\jip'^^ = Xni'l^, where ip'^^ = 
l^Cip^^. Orthogonality of ipx^j and ip'^^ may be established from the antisymmetry of 75C 
and implies that ipx^j and ip^u linearly independent eigenvectors of D\j. 

The eigenvectors ip±\u of Dfj may be expressed as linear combinations of iIj±\, with 
coefficients that depend on A and mg. When A = 0, the eigenvectors of Dh are simply given 
by ipxjj = ip\- However when A 7^ 0, they and are given by: 

i^±\H = I \. a/-zA + rUg i)x ± y^iX + nig ^-x ■ (45) 

In this basis, the matrix elements of 75 are given by: 

{>^'hH\\h) = [mg5x'^,x„ + i\\\5_x',,xu\ ■ (46) 
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Note in particular, that in the hmit A ^ rrig, the eigenstates of Dh are approximately given 
by '4'±Xh ~ V^P±'4'\ and thus they become near eigenstates of 75 with eigenvalue ±1. On 
the other hand, in the limit A ^ rUg the chirality operator is off-diagonal in this basis for 
modes which are not zero modes. 

The Hermitian Dirac operator also satisfies the commutator relation: 

[D^,75]^o = 0, (47) 

where ipo represents one of possibly many zero modes of Djj. Eq. HT] implies that 75 may 
be diagonalized within the subspace of zero modes and has eigenvalues given by ±1. Hence, 
the zero modes of and Djj are also chiral modes. On a given background gauge field 
configuration the four dimensional Dirac operator D"^ has an index given by 2Nci' = n^—n_, 
where n+ and n_ are the number of left- and right-handed chiral zero modes of D^, and the 
winding number u according to the Atiyah-Singer index theorem is given by: 

^ = 3^ / ^^'^/^"^^/^"^ ' 

where F^,^ = e^^^pF^-p. Note that since n+ — n_ is an even integer, u may be a rational 
number of the form: k/Nc, with A; G Z. If periodic BCs are used, fractional values for the 
winding number are permitted for this theory |49 ]. 

The five dimensional Hermitian DWF Dirac operator is given by Dh = TsD, where D is 
the DWF Dirac operator defined in Eq. [3l Dh satisfies a modified form of Eq. HH namely, 

{T,CyDH{T,C) = D*H . (49) 

Hence, with eigenvalues and eigenvectors oi Dh given by: 

Dh'^Kh = Ai^^An , (50) 

it similarly follows that \E'^^ = T^C^\^ are eigenvectors of Dh with eigenvalues Kh- Due to 
the antisymmetry of CFs, one may also verify that ^'a^^ and are linearly independent 
and orthogonal vectors. 

In this paper we study several quantities which may be extracted from the five dimensional 
eigenvectors of Dh- We may define a four dimensional norm: 

-Mv^l^) = $^^L(^'^)^A«(^'^) ' = 1 , (51) 
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which characterizes the profile of each eigenstate of Dh in the fifth dimension. The low 
energy modes of Dh which describe the four dimensional effective theory will appear as 
bound states, whereas at high energies, the modes will appear as propagating waves in the 
fifth dimension. We may also define a physical chirality operator Tg whose matrix elements 
are given by: 

{A'h\Ts\Ah) = J^sgn - vl/[^(x,s)vI/A,(x,s) . (52) 

x,s ^ ^ 

Near the continuum limit, matrix elements of which involve the low lying modes should 
exhibit all of the properties of 75 in the continuum four dimensional theory. 

In order to extract the residual mass and four dimensional eigenvalue density p(A), we 
consider the valence mass (m^,) dependence of the eigenvalues of Djj^ on background gauge 
field configurations, which have been generated at the input gluino mass value rrif. For 
small , we may parameterize the rriv dependence of using the reparameterized Taylor 



expansion 



48|: 



Ai,,(m,) = nl {X^ + (m, + dm,)') + 0{ml) , (53) 

and determine eigenvalue by eigenvalue the best-fit values of n^^i, Aj and 6mi, where i labels 
2-th eigenvalue of Dh ■ One may show from Eq. HQ], that with this parameterization of 
evaluated at = the gluino condensate may be expressed at order O(m^) as |:48:]: 

Note that the normalization factor 71,5 j drops out of this relation. 

We may then compare the result of Eq.[54]to the analogous continuum expression obtained 
from the operator: 

]p + mf + mres + Cs^T;^aj(^^,F;i^ + ... , (55) 

which is extracted from the Symanzik effective Lagrangian defined in Eq. [2H and where we 
have used the identity (T^dj)'"^ — f"''"^- Assuming irrelevant operators appearing in Eq. [2^ and 
consequently Eq. [55] are negligible, the form of Eq. [51] allows us to interpret the parameters 
Aj as the four dimensional eigenvalues of ip, and 5m j as an effective residual mass attributed 
to the i-th eigenvalue on a given background gauge field configuration. We may furthermore 
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estimate the eigenvalue density p'(A; nig) for the ensemble via the formula: 

where SN{X) is the number of eigenvalues within the interval A and A + (5A obtained from 
measurements on Nconf gauge field configurations. 

Close to the continuum and chiral limits, we expect for the lowest lying modes that 
the distribution of Snii values will be localized around the value for rrires, and that the 
fluctuations in 5mi around controlled by the size of finite lattice spacing artifacts. 

If on the other hand irrelevant operators in Eq. [2l] are non- negligible, the interpretation of 
the fit parameters ris j, Aj and 5mj becomes less clear. The valence Lg dependence will no 
longer reveal itself only in the parameter 5mj, but also in the parameters Aj and ns^j, due to 
the presence of the chiral symmetry breaking dimension five operators whose coefficients are 
controlled by the size of the fifth dimension. Hence, Aj obtained from Eq. [53] may no longer 
be identified as an eigenvalue Up, but rather of some more complicated operator involving 
higher order terms in the lattice spacing. 

III. THERMALIZATION AND AUTOCORRELATIONS 

For each 8^ x 8 ensemble, a total of 700 to 1000 trajectories were generated starting from 
an ordered configuration, where all gauge link fields where set to unity. For these ensem- 
bles, equilibrium was achieved within the first 250 trajectories. Approximately 1500 to 3000 
trajectories where generated for the 16^ x 32 ensembles from an ordered start and thermal- 
ization was achieved within the first 500 to 700 trajectories. Specific thermalization times 
for these ensembles are indicated in Table IIIII and Table IIVI Measurements of observables 
were made using uncorrelated configurations which were generated thereafter. A plot of the 
gluino condensate time history is shown for a variety of couplings and Lg values in Fig. [1] 
and Fig. [2] on 16^ x 32 lattices. 

In order to correctly assess the statistical errors associated with various observable mea- 
surements, we first determine the integrated correlation time Tint associated with the observ- 
able O. This may be obtained from the normalized auto-correlation function p(t), defined 
by: 

(57) 

^ iV-Ar ^ cr2 \ J 
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where 



1 ^ 
N ^ 



T 1 




(58) 



T = l 



T=l 



The integrated correlation time associated with the observable is computed using: 



,int 




(59) 



T=l 



Figure [3] and Fig. H] show the dependence of the auto-correlation function associated with 
the gluino condensate and average plaquette respectively. The exponential correlation time 
can be estimated from the value of Ar at which p(Ar) ~ e~^. The integrated correlation 
time as a function of Tmax is plotted for these same observables in Fig. [5] and Fig. [6l Be- 
cause the gluino condensate was measured using a stochastic estimator, one might expect 
decorrelation of this observable to occur relatively quickly in comparison to the average pla- 
quette. However, the tmax dependence of the integrated correlation time suggest that longer 
correlations are merely obscured by the presence of random noise. 

To better understand the correlations, we block our data with block size and consider 
the dependence of the errors on each observable. In Fig. [7] and Fig. [8] we plot the error 
as a function of block size for the gluino condensate and average plaquette. We find that for 
the gluino condensate, the decorrelation time is indeed much longer than what is implied by 
Fig. [3l Based on these considerations, we determine the gluino condensate and plaquette 
using a block size A^g = 50, where measurements have been made on every trajectory. For 
measurements of the static quark potential, we use a block size A^g = 4, where Wilson loops 
have been measured on every fifth trajectory. Unless otherwise noted, all other ensemble 
averages were performed using measurements made on every five trajectories. 

IV. MEASUREMENTS, ANALYSIS AND RESULTS 
A. Residual mass 

Measurements of the residual mass where obtained from Rif) using Eq. [35] with wall 
gluino sources. The residual mass was computed from a constant fit of -R(t) over the plateau 
region: 10-21. Fit results may be found in Table IVlT] and in Fig. [91 In the large Lg limit, 
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the residual mass is expected to behave as 

rures = ao— h — , (60) 

where the exponential term is a perturbative contribution which comes from extended modes 



above the mobility edge 



SOj. The ^ contribution to m^es niay be attributed to lattice 



dislocations and is proportional to the density of near unit modes of the five dimensional 
transfer matrix. For these simulations, the residual mass is roughly 5-20 times that of the 
input gluino mass. Figure fTOl shows a plot of the residual mass as a function of the coupling. 
The strong coupling dependence of ggests that the dislocation contribution to the 

residual mass dominates over that of the perturbative part. 

Figure [TT] shows the dependence of the residual mass on at fixed coupling and mf. 
For (3 = 2.3 and "rnj = 0.02, we fit the residual mass as a function of Ls using the functional 
form of Eq. [601 Fits were performed for a variety of different fit ranges in order to estimate 
the systematic errors associated with the fit. Results may be found in Table IVlnl These fits 
yield a negative value for the coefficient oq, which is inconsistent with the naive expectation 
that Oq should be positive. The naive expectation, however, assumes that we are simulating 
the theory at couplings which are sufficiently close to the continuum limit so that the low 
energy identity Eq. [30 holds. Provided the underlying assumptions made in the calculation 
of mj,gg(mj) are correct, one possible explanation for a negative value for oq is that higher 
order contributions to Eq. [2l|have contaminated our estimate of the residual mass obtained 
from Eq. 



B. Static potential 

Here we describe the analysis of our static potential measurements. To begin, we first 
determined the plateau region of Eq. [36] over which excited state contamination becomes 
negligible by investigating the effective potential Eq. [37|as a function of distance for a variety 
of time slice values. Figure [12], Fig. [13] and Fig. [T^show the distance dependence of \^(x, t) 
for each of the three values of and with Lg = 16 and fnf = 0.02. For /5 = 2.3, we take 
the plateau region to be the time range: 4-8, whereas for j3 = 2.353 and 2.4 we use the 
ranges: 5-9 or 5-10. For fixed distances |x|, the Wilson loops were then fit as a function 
of time to the formula Eq. [36] within the previously determined time interval. Finally the 
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extracted potentials V^(x) were fit as a function of distance to the Cornell potential given by 
Eq. [38l Table IIXI provides the fit results and respective errors which were determined by a 
jackknife analysis. Systematic errors in this fitting procedure may be estimated by varying 
the upper and lower limits of the fit ranges. These yield variations in the Sommer scale, 
however, which are comparable-if not smaller than-the statistical errors. 

Figure [15] shows a plot of the static potential at the three different values of the coupling; 
the corresponding fit parameters may be found in Table IIXI From a and a we determine the 
Sommer scale defined by Eq. [391^ As indicated in Table IIXI the Sommer scale ranges from 
between tq 3.3 for 13 = 2.3 and tq ~ 5.2 for f3 = 2.4. Our static potential results provide 
evidence that the theory is confining, which is consistent with expectations. Figure [TBI shows 
a plot of To as a function of nif + m'j.^^{mf) ^ nig for (3 = 2.3. From this plot it appears 
that the Sommer scale has little discernible gluino mass dependence at the current level of 
statistics. 



C. Gluino condensate 

In this section we describe the chiral limit extrapolation of the gluino condensate. This 
extrapolation requires taking both the Lg oo and mj ^ limits in some fashion. We 
begin by discussing the former. 

For fixed m^, we fit the gluino condensate defined in Eq. HOjas a function of Lg using the 
best available, theoretically motivated fit formula: 

ho + h——, (61) 

which is based upon the functional dependence of m^es on Lg. Note, however, that the 
dislocation contribution to m^es which appears in Eq. [60] is absent in Eq. [61] The reason 
for this omission may be understood by observing that the chiral condensate is dominated 
by contributions from UV modes, whereas the dislocation term appearing in irires niay 



predominantly be attributed to low energy phenomena 511] ■ 

Before considering the chiral limit extrapolation of the condensate we investigate the Lg 
dependence of the gluino condensate for (3 = 2.3 and a fixed input gluino mass rrif = 0.02. 



^ Measurements of the Sommer scale were performed in collaboration with I. Mihailescu. 
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The primary purpose of this investigation is two-fold. First, we would like to test our 
assumptions leading to Eq. [61] by checking that the fit value for 62 is consistent with that of 
ai obtained in Eq. EOl Second, we wish to ascertain the systematic errors associated with an 
Lg —>■ 00 extrapolation of the condensate, which come from the limited range of Lg values 
used in the fit. Results of our fits to Eq. EDare provided in Table 1X1 and Fig. [T71 From these 
fit results we conclude that the fit over the Lg range: 16-28 overestimates the Lg 00 value 
of the chiral condensate at fixed by approximately 25% compared to our most reliable 
estimates which are obtained from the fit ranges: 20-48 and 24-48. We find large variations 
in the parameter 62 with the Lg fit range, however the results are within a factor of two of 
Oi obtained from the residual mass extrapolations. 

As an alternative approach, one may consider a linear extrapolation of the gluino con- 
densate as a function of the residual mass, as was performed in js^, in order to obtain an 
Lg —* 00 value for the condensate at fixed mj. However such an approach assumes the same 
functional dependence on Lg for both the condensate and the residual mass. While indeed 
for sufficiently large Lg, each of these quantities will scale as Lj^ and such a scenario may 
be achieved, for the values of Lg used in this study, it is believed that the chiral condensate 
is dominated by a UV divergent term which is proportional to the perturbative contribution 
to the residual mass. Here, we compare the results of a linear extrapolation of the conden- 
sate as a function of the residual mass using the three data points which lie closest to the 
Lg = 00 limit (i.e., Lg = 32, 40 and 48). Results from this linear extrapolation are plotted in 
Fig. [TSl and yield a condensate value of 0.00197(4) at rrif = 0.02. Also plotted in Fig. [TSl are 
parametric plots of the gluino condensate as a function of the residual mass using Eq. [60, 
Eq. [nH and the fit parameters obtained in Table IVIIII and Table [Xl Specifically we use the 
Lg range: 24-48 fit parameters for rrires and the Lg range: 24-48 and 28-48 fit parameters for 
the condensate. A linear extrapolation of the condensate as a function of the residual mass 
yields a value which underestimates our Lg extrapolation by approximately 100-150%. 

Finally we consider the chiral limit extrapolation of the condensate. Such an extrapolation 
may be achieved by considering the individual Lg and rrif dependence of the condensate, as 



was considered in 



31| . Here, we perform chiral limit extrapolations of the gluino condensate 



at a single lattice spacing {(3 = 2.3) using two different limit orders. First, we perform a 



23 



linear m/ 



extrapolation of the gluino condensate at fixed Lg using the formula 

Co + Cim/ , 



(62) 



followed by an ^ oo extrapolation of the rrif = result using Eq. [6T] We believe that this 
approach is more transparent than, for instance, a linear extrapolation of the condensate as 
a function of rrires (followed by a linear — > extrapolation), and most importantly, it 
correctly accounts for the functional form of Eq. [6T1 

Following the double extrapolation procedure outlined above, we obtain an unrenormal- 
ized value for the gluino condensate in the chiral limit at finite lattice spacing. Results from 
these fits are tabulated in Table |XT] and Table Kill and plotted in Fig. [19] and Fig. [20|, and 
yield the chiral limit value of 0.00320(9) for the condensate. Chiral extrapolations have also 
been performed by reversing the order of limits and yield consistent results with comparable 
statistical error bars. The fits of the condensate as a function of Lg were performed over a fit 
range: 16-28 and, as such, we expect our chiral limit result to over estimate the true answer 
by approximately 25%, as was found for our extended Lg fits at /3 = 2.3 and rrif = 0.02, 
which are described above. We have performed additional fits using other, phenomenolog- 
ically motivated fit formulae to describe the Lg dependence of the gluino condensate (e.g., 
Eq. EH without the Lj^ prefactor in the second term). Such fits yield an approximate 20% 
variation in the chirally extrapolated value of the gluino condensate as compared to that 
obtained by using Eq. [6TJ 

D. Spectrum of the Hermitian DWF Dirac operator 

Here we study the eigenvalues and eigenvectors of the Hermitian DWF Dirac operator. 
The primary purpose of this study is to establish to what degree the five dimensional theory 
reproduces the correct four-dimensional, continuum low energy physics and to obtain an 
independent measure of the residual mass and gluino condensate. To achieve this aim, we 
measure the lowest 64 eigenvalues and eigenvectors of the Hermitian DWF Dirac operator 
on 8"^ X 8 lattices using the method of Kalkreuter-Simma (KS) [52]. In this method, the 
Hermitian matrix Djj is first diagonalized using a modified Rayleigh-Ritz diagonalization 
procedure, where we have exploited the relationship between degenerate eigenvectors \1/a^ 
and in order to eliminate any unnecessary minimizations of the Ritz functional. Subse- 
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quently, the matrix Dh is diagonalized using Jacobi's method within the subspace obtained 
from diagonahzing Djj. These steps are then iterated until specified stopping criteria have 
been achieved. A detailed description of the procedure used to diagonalize and the stop- 



ping criteria may be found in 53|. For our purposes, the process was considered converged 
when the change in the eigenvalues between iterations was less than 1 x 10^^. We check 
this choice by increasing the stopping condition to 1 x 10"^'' on a test configuration and 
determined that eigenvalues changed by at most 0.02% over the the range of eigenvalues 
considered. 

The eigenvectors of Dh obtained from the KS method are ordered such that the magni- 
tude of Ah are ascending with eigenvector number. In all of the analysis that follows we use 



only the first 60 of the 64 eigenvalues obtained, since the last few may be unreliable |53 |. 
Since there is a two-fold degeneracy in the spectrum due to Eq. HU this yields a total of 30 
non-degenerate eigenvectors and eigenvalues. All measurements in this work were performed 
on 8^ X 8 lattices with f3 = 2.3 and a sea gluino mass: trif = 0.02. In physical units, the 
lowest 60 eigenvalues for a typical configuration range roughly between < ro|A£f| < 0.07 
for the ensembles and valence masses under consideration. 

The lowest 60 eigenvectors were used to compute the matrix elements of the physical 
gamma matrix defined in Eq. [52l as well as to determine the localization of wave functions 
in the fifth dimension, which may be characterized by A/'a^(s) as defined in Eq. [511 Figure [2T] 
and Fig. [22] shows the results of these calculations for a representative background gauge field 
configuration in the ensemble: f3 = 2.3, Ls = 24 and mj = 0.02, with m„ = 0.02. As may be 
seen in these figures, the eigenvectors are exponentially localized on both boundaries of the 
fifth dimension and alternate with increasing Kh- The near-diagonal nature of F^ confirms 
that the low-lying eigenvectors of Dh are near chiral modes, however these states are not 
necessarily near-zero modes of Dh- We find that in every instance: 

sgn(AH) =sgn((AH|F,|AH)) , (63) 

which is entirely consistent with continuum predictions in the regime where A < m^, accord- 
ing to Eq. [46l 

Next we study the my dependence of A^ and attempt to extract the residual mass and 
chiral condensate for the /3 = 2.3, Ls = 16, 20 and 24, and frif = 0.02 ensembles. Eigenvalues 
were computed for five values of m„, ranging from —0.16 to —0.12, on a total of 150 gauge 
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configurations. After appropriately reordering eigenvalues in order to account for possible 
level crossings, we fit the lowest 60 eigenvalues of Djj^ as a function of to Eq. [53] and 
extract the parameters ns j, Aj and Snii for each eigenvector. Figure [251 shows a plot of the 
fit results for the lowest 20 eigenvalues on a representative background field configuration 
for Lg = 24. As may be seen in this plot, for a typical eigenvector, the splitting of iA// 
pairs due to the finite lattice spacing effects appear comparable to the spacings between 
adjacent eigenvectors. As a qualitative measure, we conclude that the lattice spacing is not 
sufficiently small for the emergence of this continuum behavior of the Dirac spectrum. 

In Fig.[24]we plot the distribution of 6mi values at fixed eigenvalue number for the Ls = 24 
ensemble. Identifying the peak of the distribution of 6mi values with the residual mass relies 
on an assumption that the spectrum of Dh is in some sense continuum-like. Unfortunately 
due to the lack of eigenvalue pairing, such interpretation is somewhat obscured. For the 
lowest eigenvalue, the mean of the distribution appears consistent with the value of rrires 
obtained from Eq. [HH however for larger eigenvalues, the mean increases in magnitude. 

In Fig. [25] we plot the values of A, verses 5mj for all 60 x 150 eigenvalues, which were 
extracted from fits to Eq. [53] for the Lg = 24 ensemble. At low energies, it is evident 
that we are in a regime where Aj < Srrii, which is consistent with our findings for the 
structure of in Fig. [22] If the distribution of 6mi is highly peaked about nires, one may 
in principle tune ~ —rrires such that rrig < A and expect off diagonal pairing to emerge 
in the matrix elements of F^. An attempt to produce off-diagonal pairing of Fs by tuning 
m„ ^ —m'r^gimf) was unsuccessful, but is likely due to the fact that off-diagonal pairing in 
Vg may only occur when the ± eigenvalues of Dh are paired. In Fig. [26] we plot Nk^^s) 
for the same configuration as in Fig. [JT] however with = —0.25. Figure [27] displays 
the corresponding matrix elements of Fg for this configuration. From the first plot, we see 
that the lowest lying wave functions are no longer strongly localized on the fifth dimension 
boundaries. As such, the diagonal structure of F^ in Fig. [27] disappears for the lowest few 
modes, although off-diagonal pairing in not very evident either. 

In Fig. [28] we plot the density of four dimensional eigenvalues p'{X,mg) for Lg = 16, 20 
and 24, which have been extracted from 60 x 150 eigenvalue measurements using Eq. [53] 
The distribution exhibits a relatively constant region between A ~ 0.15 and 0.25, and a 
depleted region between A ~ and 0.15 which may be attributed to finite volume effects. 
Contrary to continuum expectations, there is no visible peak at A = 0, which one would 
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normally attribute to near-zero modes. If such a peak were to exist, it may be that its 
height is suppressed due to a finite bin size 6X, which is too large. From Fig. [231 however, it 
appears that the number of zero modes is few, suggesting that there is perhaps little topology 
change and that we are primarily within the zero topological charge sector of the theory. 
This observation is consistent with the findings of [31], where fiuctuations in the condensate 
as a function of MD time where studied as a possible indicator of topology change. In 
that study, no large spikes in evolution of the condensate (which would correspond to the 
presence of zero-modes) were observed for the same choice of parameters. Yet in light of the 
fact that the residual mass is so large in these simulations, it may also be that the fermions 
simply provide a poor measure of topology. To make a more definitive statement regarding 
the topic, it would be interesting and beneficial to study the gauge fields directly instead. 

In principle one may extract an estimate of Pgjj(0,mg) from the constant region of the 
eigenvalue density, compute the condensate and finally perform an nig extrapolation. 
But for the current choice of simulation parameters, the low lying spectrum appears to be 
too heavily distorted by the large residual mass for such an analysis to be reliable. We may 
still extract some useful information from the eigenvalue density, however. In our analysis 
of the chiral extrapolation of the gluino condensate described in Sec. IIV C\ we assumed that 
the Ls dependent contribution to the condensate was dominated by UV modes. We may 
check this assumption, by computing the contribution to Eq. HQ] from modes in the interval: 
0-A, using our fit values for Aj and Srrii and Eq. [54] with = rrif. In Fig. [29] we plot the 
integrated condensate as a function of the upper limit value A for Ls = 16, 20 and 24. The 
contribution to the condensate from modes |A|ro < 0.8 is approximately 1/3 of the total 
value for each choice of Lg, which is consistent with our assumptions. 

Finally for completeness. Fig. [30] provides a histogram of values extracted from the 
A^(m^) fits for the (3 = 2.3, Lg = 24 and mj = 0.02 ensemble, where it is evident that 
is strongly peaked at approximately ~ 0.16. 

V. CONCLUSION 

We have performed dynamical numerical simulations of N=l SYM theory on 8^ x 8 and 
16^ X 32 lattices using domain wall fermions with several goals in mind: to establish a 
lattice scale, assess the size of chiral symmetry breaking effects at finite Lg, and understand 
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the degree to which DWFs are correctly reproducing the desired continuum physics at low 
energy. In part, these efforts were motivated by the necessity to establish benchmarks for 
future studies. As such, our work consists of a variety of basic measurements on 16^ x 32 
lattices, including measurements of the static potential, residual mass and gluino condensate. 
We buttress these results by analyzing the eigenvalues and eigenvectors of the Hermitian 
DWF Dirac operator on 8^ x 8 lattices, from which an independent estimate of the residual 
mass and chiral condensate may in principle be extracted. We briefly summarize each of 
our measurements and results below. 

We have presented results for the residual mass, which have been obtained using two 
different methods. The first approach utilized a low energy identity relating the wall pseudo- 
scalar density to the midpoint pseudo-scalar density and assumptions about the decompo- 
sition of connected and disconnect parts of various pseudo-scalar correlators. Using this 
approach we obtained residual masses which where approximately 5-20 times larger than 
the input gluino mass m/. A second approach for obtaining m^es, based upon the rriy de- 
pendence of the eigenvalues of D'jj, gave inconclusive yet qualitatively consistent results for 
the residual mass. On an eigenvalue by eigenvalue basis, chiral symmetry breaking shifts in 
nif where comparable, if not larger than the value for m^es obtained in the first approach. 
The limited success of the second method in determining a precise value for the residual 
mass may be attributed to large finite lattice spacing artifacts. 

The static potential for fundamental representation sources was determined using stan- 
dard techniques for three different values of the couphng. The potential exhibits characteris- 
tics indicative of confinement. Prom the static potential we obtain estimates of the Sommer 
scale, which range from approximately tq ^ 3.3 at /9 = 2.3 to tq ^ 5.3 aX f3 = 2.4. 

We study properties of the low energy eigenvalues of the Hermitian DWF Dirac operator 
in an effort to understand to what degree the desired low energy physics is reproduced by 
the DWF formahsm. By studying the profile of the wave functions as a function of the fifth 
dimension coordinate, we have established that the low energy modes are indeed localized on 
the right and left boundaries of the fifth direction. And, the matrix elements of the physical 
chirality operator Tg arc consistent with the presence of a large gluino mass compared to 
the typical low energy eigenvalues of ip. 

We have performed a chiral extrapolation of the gluino condensate at a fixed value of 
the lattice spacing, /3 = 2.3. The chiral extrapolation of the condensate was performed in a 
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variety of different ways in order to elucidate some of the systematic errors inherent with the 
extrapolation. For the parameters under consideration we find that various extrapolation 
procedures yield, at best, a value of the condensate which is reliable to approximately 25%. 
We attempt to provide an independent measurement of the gluino condensate from our 
studies of the Dirac spectrum. However, given the current large values of the residual mass, 
a reliable estimate of the gluino condensate from the Banks-Casher relation was not possible. 

Finally, we come to the general conclusion that there is, at present, no evidence to suggest 
that our simulations are not in the same universality class as SYM at non-zero gluino mass. 
However, additional simulations will be required at far smaller residual masses and weaker 
couplings in order to achieve a more reliable measure of the quantities discussed in this paper. 
Reducing the residual mass may be achieved with a variety of approaches; the simplest 
approach is to increase the size of the fifth dimension. This method for reducing the residual 
mass becomes costly, however, once the dislocation term appearing in Eq. [60] dominates over 
the perturbative contribution to the residual mass. Alternatives to this approach, which may 
be more efficient, include going to weaker coupling, or using an improved action such as the 
DBW2 5J] , Iwasaki 55|, l56[|, o r auxiliary determinant 57l| a ction. The latter action combines 



the Gap DWF method of 
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59| with the approach of 



601 1 in order to suppress the residual 



mass while still maintaining adequate topological tunneling. Such considerations will be left 
to a future study. Finally we emphasis that scaling analysis at weaker coupling is essential in 
order to truly understanding the effects of finite lattice spacing on observables, and whether 
or not the gluino condensate survives the continuum limit. 
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APPENDIX A: CORRELATORS INVOLVING MAJORANA FERMIONS 

In A/" = 1 SYM, the fundamental fermionic degrees of freedom are Majorana fermions 
rather than Dirac fermions. As a resuh, correlation functions written in terms of propa- 
gator contractions may differ from that of conventional QCD. Here we briefly review the 
construction of some of the basic observables studied in the paper, although for simplicity 
we limit the majority of the discussion to four dimensions. Generalization of the results to 
the specific case of DWFs will be considered at the end of this section. 

To begin with, consider the generic Majorana fermion path integral: 



where the indices for k = 1, . . . , 2N collectively represent space-time coordinate (x), spin 
(a) and color (c), and summation over these indices is implied. Due to the anti-commuting 
nature of the variables ip, components of Gij^...i2N which are symmetric under the interchange 
of any two indices will not contribute to G^^^^ The Majorana matrix M is related to the 
Dirac operator via the relation M = CD^, where C is the charge conjugation matrix 
whose properties in four dimension are given in Eq. [9l We allow the Dirac operator and Q 
to depend on background field configurations and, as such, so may the normalization factor 
A. By definition we choose A such that C-^^ = 1, and so it follows that A = Pf(M). 
Expectation values of operators involving the fermion fields ip in the full gauge theory are 
given by gauge averages over G^"^^^ with the probability measure A e~^° , where Sq is the 
gauge action. 

The path integral given by Eq. lAll may be evaluated in a standard fashion, and yields 




(Al) 



'11...12N ' 



(A2) 



where 



m; 



r-l 



-1 \ 



Pf 



(A3) 




M. 



30 



is the Pfaffian constructed from Majorana propagators = [D^]~^C\ which start and 
end on the points ii . . . i2N in ah possible combinations. Note that since the matrix M is 
antisymmetric, the diagonal matrix elements of the propagator vanish. 

As an example of how this result may be applied, consider the two important cases = 1 
and iV = 2, for which: 



= Qi^^iiMr} Mr} - Mr} Mr} +Mr}Mr}) . (A4) 



In the case of the condensate — (V'V') defined in Eq. HQl we have: 

which according to the expression for ^(2) in Eq. El yields: 

^5^Tr^(x,x) (A6) 

X 

prior to gauge averaging. Here, S = [D^]~^ is the Dirac propagator and the trace is taken 
over both spin and color. In the case of the pseudo-scalar correlator (P(x)P(x')), where 
P{x) = iIj{x)'J^i/j{x), we have 

^ii«2i3*4 ~ ^xi,x^X2,x^X3,x'^X4,x' i.C'y^') ai,a2i.^'~f5) a3,a3^ci,C2^a3,a3 ! (-'^'^) 

which according to the expression for G(^) in Eq. |Al yields: 

2Tr [S'(x,x')755'(x',x)75] - Tr [S'(x,x')75] Tr [S'(x',x)75] , (A8) 

prior to gauge averaging. Note that the connected and disconnected parts of this correlator, 
as defined in Eq. [321 correspond to the single and double trace terms respectively. 

While the results presented here pertain to the four dimensional theory, they may easily be 
generalized to the case of DWFs, where the interpolating gluino fields q{x) are expressed in 
terms of the fields s) at the fifth dimension boundary, according to Eq. [TTl Specifically 
the results for the gluino condensate and pseudo-scalar correlator given by Eq. IA6I and 
Eq. I A8 1 remain valid with the simple replacement S — > Sq, where 

Sq{x,x') = P-D^^Q.y j^^_^P^ + P_P)^ J.^, oP+ + 
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and D is the DWF Dirac operator defined in Eq. [31 
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TABLE I: Parameter values for the RHMC algorithm for 8^ x 8 lattices. Xmax and Xmin are the 
maximum and minimum eigenvalues of V^M^jUif) required for the rational approximation, and 
umd and umc represent the degree of the rational approximation for the MD and MC accept/reject 
step, respectively. 



/? Ls nif 


Xmax 


Xmin 




riMD nMC 


\PV 
^max 


\PV 
min 




„PV 
"'MD 


„PV 


2.3 12 0.02 


4.0 


4 X 10" 


-4 


9 


15 


4.0 


1 X 10- 


-3 


6 


9 


16 0.02 


4.0 


4 X 10" 


-4 


9 


15 


4.0 


1 X 10" 


-3 


6 


9 


20 0.02 


4.0 


4 X 10" 


-4 


9 


15 


4.0 


1 X 10" 


-3 


6 


9 


24 0.02 


4.0 


4 X 10" 


-4 


9 


15 


4.0 


1 X 10" 


-3 


6 


9 
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TABLE II: Parameter values for the RHMC algorithm for 16^ x 32 lattices. An explanation of 
these parameter values may be found in Table HI 



/3 


Ls ruf 


^max 




A 


min 






nMC 


\PV 

^max 




\PV 

vain 




"'MD 


"'MC 


2.3 


16 0.01 


4.0 


4 


X 


10- 


-4 


9 


15 


4 





1 


X 


10- 


-2 


6 


9 




0.02 


4.0 


4 


X 


w 


-4 


9 


15 


4 





1 


X 


10- 


-2 


6 


9 




0.04 


4.0 


4 


X 


10- 


-4 


9 


15 


4 





1 


X 


10- 


-2 


6 


9 




20 0.01 


4.0 


4 


X 
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15 


4 
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X 


10- 


-2 


6 


9 




0.02 


4.0 


4 


X 
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-4 
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15 


4 





1 


X 


10- 


-2 


6 


9 




0.04 


4.0 


4 


X 


10" 


-4 


9 


15 


4 





1 


X 


10- 


-2 


6 


9 




24 0.01 


4.0 


4 


X 
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-4 
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15 


4 





5 
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10- 


-3 


6 


9 




0.02 


4.0 


4 
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-4 
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15 


4 





5 


X 


10- 


-3 


6 


9 




0.04 


4.0 


4 


X 


10- 


-4 


9 


15 


4 





5 


X 


10- 


-3 


6 


9 




28 0.01 


4.0 


4 
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-4 
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15 


4 





5 


X 


10- 


-3 


6 


9 




0.02 


4.0 


4 


X 
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-4 


9 


15 


4 





5 


X 


10- 


-3 


6 


9 




0.04 


4.0 


4 


X 


10- 


-4 


9 


15 


4 





5 


X 


10- 


-3 


6 


9 




32 0.02 


4.0 


1 


X 


10- 


-4 


9 


15 


4 





1 


X 


10- 


-3 


6 


9 




40 0.02 


4.0 


4 
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10- 


-5 
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4 





5 
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10- 


-4 


6 


9 




48 0.02 
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X 
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-5 


9 


15 


4 





5 


X 


10- 


-4 


6 


9 


2.353 


16 0.02 


4.0 


4 X 10" 


-4 


9 


15 


4.0 


1 X 10- 


-2 


6 


9 




28 0.02 


4.0 


4 X 10" 


-4 


9 


15 


4.0 


5 X 10- 


-3 


6 


9 


2.4 


16 0.02 


4.0 


4 X 10" 


-4 
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15 


4.0 


1 X 10- 


-2 


6 


9 




28 0.02 
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4 X 10" 


-4 
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15 


4.0 


5 X 10- 


-3 
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9 
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TABLE III: MD evolution parameters and MC statistics for 8^ x 8 lattices. For each ensemble, the 
MD evolution trajectory length was 5x St MD time units. A total number of trajectories for each 
ensemble is given by Ntraj- Nthf,rm is the number of initial trajectories before configurations in a 

given ensemble were deemed thermalizcd. "Accept" indicates the acceptance rate associated with 
the MC accept /reject step. Ensemble averages of various functions of the Hamiltonian H were 
performed using m(\xsmxnii(nil s on (nTry Irajocl.ory. 



f3 Ls nif 


5r 


Ntraj 


^therm 


Accept {AH) y/{AH^) (e"^^) 


2.3 12 0.02 


0.260 


1000 


300 


0.626(18) 0.498(28) 1.095(27) 0.956(39) 


16 0.02 


0.220 


750 


300 


0.749(20) 0.146(26) 0.577(19) 1.004(19) 


20 0.02 


0.240 


1000 


300 


0.594(18) 0.541(39) 1.165(30) 1.013(60) 


24 0.02 


0.230 


1000 


300 


0.627(18) 0.470(35) 1.965(43) 0.966(42) 



37 



TABLE IV: MD evolution parameters and MC statistics for 16^ x 32 lattices. An explanation of 
these parameter values may be found in Table IIIII 



13 Ls nif 






^therm 


Accept 


{AH) 






2.3 16 0.01 


0.160 


3125 


500 


0.758(19) 


0.193(28) 


0.673(9) 


1.016(14) 


0.02 


0.160 


3195 


500 


0.762(18) 


0.187(26) 


0.625(8) 


0.991(12) 


0.04 


0.160 


2790 


500 


0.776(19) 


0.152(26) 


0.577(8) 


1.003(12) 


20 0.01 


0.155 


2895 


500 


0.745(20) 


0.208(30) 


0.685(10) 


1.008(15) 


0.02 


0.155 


2655 


500 


0.753(21) 


0.212(31) 


0.674(10) 


0.993(16) 


0.04 


0.160 


2760 


500 


0.731(21) 


0.238(32) 


0.722(11) 


0.995(16) 


24 0.01 


0.155 


2855 


500 


0.775(19) 


0.142(26) 


0.578(8) 


1.019(13) 


0.02 


0.130 


2620 


500 


0.792(20) 


0.143(28) 


0.593(11) 


1.020(14) 


0.04 


0.155 


2610 


500 


0.760(21) 


0.175(32) 


0.690(10) 


1.050(17) 


28 0.01 


0.135 


2740 


500 


0.817(18) 


0.086(22) 


0.474(7) 


1.022(11) 


0.02 


0.140 


2855 


500 


0.796(19) 


0.158(24) 


0.536(7) 


0.974(11) 


0.04 


0.155 


2880 


500 


0.784(19) 


0.157(25) 


0.577(8) 


0.996(12) 


32 0.02 


0.140 


2730 


500 


0.757(20) 


0.180(29) 


0.654(9) 


1.016(15) 


40 0.02 


0.140 


2189 


500 


0.628(27) 


0.455(51) 


1.028(18) 


0.974(27) 


48 0.02 


0.140 


1555 


500 


0.608(34) 


0.586(76) 


1.247(25) 


1.031(46) 


2.353 16 0.02 


0.163 


2575 


600 


0.782(20) 


0.160(27) 


0.583(9) 


0.996(13) 


28 0.02 


0.140 


2305 


600 


0.829(20) 


0.072(22) 


0.415(7) 


1.012(10) 


2.4 16 0.02 


0.160 


2710 


750 


0.824(18) 


0.093(20) 


0.432(7) 


0.996(10) 


28 0.02 


0.140 


2285 


750 


0.882(18) 


0.040(18) 


0.320(6) 


1.010(8) 
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TABLE V: Gluino condensate {qq), {q'Jdq) and average plaquette (P) for 8^ x 8 lattices. 



P Ls ruf 


(qq) 


{qj^q) X 10- 


-5 (P) 


2.3 12 0.02 


0.010631(40) 


-1(6) 


0.73496(36) 


16 0.02 


0.008736(42) 


-5(9) 


0.73296(51) 


20 0.02 


0.007522(37) 


-10(6) 


0.73276(30) 


24 0.02 


0.006806(34) 


-1(7) 


0.73202(32) 
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TABLE VI: Gluino condensate {qq), {q^y^q) and average plaquette (P) for 16^ x 32 lattices. 



(3 Ls nif 


{-qq) 


{qi5q) X 10" 


-6 (P) 


2.3 16 0.01 


0.0078277(43) 


11(6) 


0.733463(33) 


0.02 


0.0086506(30) 


3(5) 


0.733218(24) 


0.04 


0.0102120(41) 


-1(5) 


0.733203(36) 


20 0.01 


0.0067032(46) 


-1(6) 


0.732618(36) 


0.02 


0.0075158(46) 


5(6) 


0.732530(57) 


0.04 


0.0091267(35) 


1(6) 


0.732419(31) 


24 0.01 


0.0059855(40) 


2(7) 


0.732127(35) 


0.02 


0.0068081(42) 


-2(7) 


0.732110(41) 


0.04 


0.0084396(36) 


1(6) 


0.731983(37) 


28 0.01 


0.0055027(46) 


-8(7) 


0.731734(34) 


0.02 


0.0063346(33) 


6(6) 


0.731688(32) 


0.04 


0.0079882(32) 


-3(5) 


0.731628(38) 


32 0.02 


0.0059947(47) 


11(7) 


0.731418(45) 


40 0.02 


0.0055266(43) 


-2(10) 


0.730947(44) 


48 0.02 


0.0052238(43) 


-10(10) 


0.730727(50) 


2.353 16 0.02 


0.0077917(72) 


7(6) 


0.745243(40) 


28 0.02 


0.0057106(61) 


7(10) 


0.743853(38) 


2.4 16 0.02 


0.0068851(64) 


-8(8) 


0.754149(30) 


28 0.02 


0.0049179(82) 


-16(9) 


0.752951(27) 
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TABLE VII: Fit results for the residual mass which was extracted from R{t) for 16^ x 32 lattices. 
The column labeled "time range" indicates the plateau region over which the R{t) was fit. 



f3 


Ls rrif 


time ran^ 


re xVd-o.f. 




2.3 


16 0.01 


10-21 


63.8/11 


0.18682(13) 




0.02 


10-21 


15.0/11 


0.18939(12) 




0.04 


10-21 


17.2/11 


0.19223(13) 




20 0.01 


10-21 


10.8/11 


0.17151(14) 




0.02 


10-21 


35.1/11 


0.17340(14) 




0.04 


10-21 


49.4/11 


0.17680(13) 




24 0.01 


10-21 


13.4/11 


0.15807(13) 




0.02 


10-21 


18.4/11 


0.15906(15) 




0.04 


10-21 


19.8/11 


0.16367(13) 




28 0.01 


10-21 


28.6/11 


0.14685(13) 




0.02 


10-21 


25.5/11 


0.14834(13) 




0.04 


10-21 


111.7/11 0.15283(13) 




32 0.02 


10-21 


5.1/11 


0.13905(14) 




40 0.02 


10-21 


23.2/11 


0.12397(14) 




48 0.02 


10-21 


29.1/11 


0.11225(16) 


2.353 16 0.02 


10-21 


45.5/11 


0.14071(19) 




28 0.02 


10-21 


42.5/11 


0.10269(18) 


2.4 


16 0.02 


10-21 


9.5/11 


0.10082(15) 




28 0.02 


10-21 


24.6/11 


0.06513(17) 
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TABLE VIII: Fit results for the residual mass as a function of Ls for 16^ x 32 lattices with /3 = 2.3 
and rnf = 0.02. The column labeled by "L^ range" indicates the range of Lg values used in the fit. 
The parameters oq, ai and 02 are defined in Eg. [5U1 



13 Lg range ruf 


xVd.o.f. 


ao ai 


^2 


2.3 16-48 0.02 


53.5/4 


-6.20(3) 0.0275(4) 


7.03(5) 


20-48 


5.9/3 


-6.37(7) 0.0236(7) 


7.43(10) 


24-48 


2.2/2 


-6.30(5) 0.0255(11) 


7.24(12) 
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TABLE IX: Static quark potential fit parameters for a subset of 16^ x 32 lattices. The column 
labeled by "time range" indicates the time window over which Wilson loops were fit at fixed 
distances in order to extract the potential. The column labeled by "distance range" indicates the 
distance window over which the potential was fit. 






J-^S ' ' ^ J 


time ran^ 


5e distance ran 


Vn Of 


cr 


Tn 
' u 


2.3 


16 0.02 


4-8 


\/2-6 


511fl0) 176fl0) 


0.132(2) 


3 344f21) 




0.02 


4-8 


\/3-\/29 


483(17) 140('23') 


0.137(4) 


3 319f22') 




0.02 


4-8 


^/3-6 


501fl8) 161f23'l 


0.134(4) 


3 339f23') 




0.04 


4-8 


a/2-6 


524fl5l 185fl4) 


129f3) 


3 376 f 29) 




0.04 


4-8 


V3-\/29 


0.516(26) 0.178(33) 


0.130(5) 


3.381(32) 




0.04 


4-8 


\/3-6 


534f26) 200f32) 


0.127(5) 


3 359f311 




20 0.02 


4-8 


a/2-6 


0.511(13) 0.174(13) 


0.134(4) 


3.322(47) 




0.02 


4-8 


a/3- a/29 


480 f44) 132 f 54) 


139fll) 


3 306 f 77) 




0.02 


4-8 


a/3-6 


483f411 136f49) 


138fl0) 


3 309f71) 




0.04 


4-8 


\/2-6 


532fl5) 197fl5) 


0.131(4) 


3 327f31) 








\/^-\/9Q 




U. ±^Ut (J I 


o.^y X \^-^ ) 




0.04 


4-8 


a/3-6 


0.507(27) 0.163(32) 


0.135(5) 


3.316(34) 




24 0.02 


4-8 


a/2-6 


0.522(19) 0.182(18) 


0.132(4) 


3.340(38) 




0.02 


4-8 


\/^-\/9Q 


n 504f39'l f) 1 ^RfSQ"! 


1 '^'^(7) 

\J. J. (JO \ 1 J 


3 33nf41 ) 




0.02 


4-8 


\/3-6 


0.511(31) 0.166(38) 


0.134(6) 


3.322(41) 




0.04 


4-8 


\/2-6 


0.524(15) 0.185(14) 


0.129(3) 


3.376(29) 




0.04 


4-8 


a/3- a/29 


0.516(26) 0.178(33) 


0.130(5) 


3.359(31) 




0.04 


4-8 


\/3-6 


0.534(26) 0.200(33) 


0.127(5) 


3.381(32) 


2.353 16 0.02 


5-9 


a/2-6 


0.543(13) 0.202(13) 


0.077(3) 


4.324(64) 




0.02 


5-9 




0.550(23) 0.218(32) 


0.077(4) 


4.313(76) 




0.02 


5-9 


V3-6 


0.569(23) 0.240(23) 


0.074(4) 


4.379(80) 


2.4 


16 0.02 


5-9 


V2-6 


0.528(7) 0.191(7) 


0.053(2) 


5.254(67) 




0.02 


5-9 


V3-V29 


0.534(16) 0.199(22) 


0.052(3) 


5.272(101) 




0.02 


5-9 


V3-6 


0.536(15) 0.202(20) 


0.052(3) 


5.288(94) 




0.02 


5-10 


a/2-6 


0.533(5) 0.196(5) 


0.052(1) 


5.280(45) 




0.02 


5-10 


V3-V29 


0.532(12) 0.196(16) 


0.053(2) 


5.257(71) 




0.02 


5-10 


V3-6 


430.539(11) 0.205(15) 


0.051(2) 


5.306(68) 



TABLE X: Fit results for gluino condensate as a function of Ls for 16^ x 32 lattices with (5 = 2.3 
and rrif = 0.02. The column labeled by "L^ range" indicates the range of Lg values used in the fit. 
The parameters 6o) and 62 are defined in Eq. [HTJ 



j3 Ls range mj 


xVd.o.f. 


bo hi 62 


2.3 16-28 0.02 


0.1/1 


0.004994(53) 0.1057(18) 0.0370(20) 


16-32 


1.9/2 


0.004939(39) 0.1040(14) 0.0350(15) 


16-40 


26.8/3 


0.004758(20) 0.0989(7) 0.0289(7) 


16-48 


62.9/4 


0.004650(16) 0.0962(5) 0.0255(5) 


20-48 


14.8/3 


0.004529(29) 0.0897(9) 0.0204(10) 


24-48 


1.9/2 


0.004392(62) 0.0839(18) 0.0154(19) 



44 



TABLE XI: Fit results for gluino condensate as a function oimj for 16^ x 32 lattices with (3 = 2.3. 
The column labeled by "mj range" indicates the range of values used in the fit. The parameters 
Co and ci are defined in Eq. 



(3 Ls rrif range 


xVd.o.f. 


Co Ci 


2.3 16 0.01-0.04 


40.6/1 


0.0070544(51) 0.07915(19) 


20 


0.7/1 


0.0058979(55) 0.08073(19) 


24 


0.8/1 


0.0051697(49) 0.08177(18) 


28 


0.5/1 


0.0046770(51) 0.08279(18) 
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TABLE XII: Fit results for the mj = extrapolated value of the gluino condensate as a function 
of Lg for 16^ X 32 lattices with (3 = 2.3. The column labeled by "L^ range" indicates the range of 
Lg values used in the fit. The parameters Bq, bi and 62 are defined in Eq. ETJ 



f3 Lg range nif 


xVd-o.f. 60 hi 62 


2.3 16-28 


0.07/1 0.00320(9) 0.1051(24) 0.0334(27) 
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Trajectory number 



FIG. 1: Gluino condensate as a function of trajectory number for 16^ x 32 lattices with /3 = 2.3 
and nif = 0.02. 
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FIG. 2: Gluino condensate as a function of trajectory number for 16^ x 32 lattices with Lg 
and nif = 0.02. 
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FIG. 3: Autocorrelation function associated with the gluino condensate for 16^ x 32 lattices with 
P = 2.3 and m/ = 0.02. The dashed line indicates the location of e~^. 
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FIG. 5: Integrated correlation time associated with the gluino condensate for 16^ x 32 lattices with 
P = 2.3 and m/ = 0.02. 
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FIG. 6: Integrated correlation time associated with the average plaquette for 16^ x 32 lattices with 
j3 = 2.3 and m/ = 0.02. 
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Block size 



FIG. 7: Gluino condensate error as a function of block size for 16^ x 32 lattices with (3 = 2.3 and 
ruf = 0.02. 
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FIG. 9: R as a function of time for 16^ x 32 lattices with P = 2.3 and ruf = 0.02. Solid lines 
indicate the value of m'^^^^mf) obtained from a constant fit to data, as given in Table rvTIl 
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FIG. 11: Residual mass as a function of Ls for 16^ x 32 lattices with (3 = 2.3 and nif = 0.02. Solid 
curves represent fit results obtained for a variety of different Lg ranges. Details of the fits may be 
found in Table IVnH 
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FIG. 12: Effective static potential Vf,ff{x,t) as a function of distance |x| at fixed time t for a 
16^ X 32 lattice witli /3 = 2.3, = 16 and m/ = 0.02. 
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FIG. 13: Effective static potential Vf,ff{x,t) as a function of distance |x| at fixed time t for a 
16^ X 32 lattice witli /3 = 2.353, = 16 and m/ = 0.02. 
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FIG. 14: Effective static potential Vf,ff{x,t) as a function of distance |x| at fixed time t for a 
16^ X 32 lattice witli /3 = 2.4, = 16 and m/ = 0.02. 
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FIG. 15: Static potential as a function of distance for 16^ x 32 lattices with Lg = 16 and ruf = 0.02. 
The solid curves represent fits to the potential given in Table ITXl and dashed lines indicate la error 
bars for the corresponding values for the Sommer scale obtained from Eq. [39j 
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FIG. 16: Sommer scale as a function of mj + m',j.^g{'m f) « rUg for 16^ x 32 lattices with /3 = 2.3. 
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FIG. 17: Gluino condensate as a function of Ls for 16^ x 32 lattices with /3 = 2.3 and rnj = 0.02. 
Solid curves represent fit results obtained for a variety of different Lg ranges. Details of the fits 
may be found in Table [Xj 
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FIG. 18: Gluino condensate as a function of the residual mass for 16^ X 32 lattices with /? = 2.3 
and mf = 0.02. The two curve labeled "parametric" represent curves obtained from the fit results 
obtained in Table IVnll and Table 1X1 The fit parameters used for the residual mass are the same, 
and given by the Lg range: 24-48. Parameters used for the condensate are indicated on the plot. 
The curve labeled "linear" represents a linear extrapolation of the Lg = 32, 40 and 48 results. 
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FIG. 19: Gluino condensate as a function of for 16^ x 32 lattices with (3 = 2.3. Solid curves 
represent fit results which are obtained from Table IXH 
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FIG. 20: ruf = extrapolated values of the gluino condensate as a function of Lg for 16^ x 32 
lattices with f3 = 2.3. The solid curve represents the fit result obtained from Table IXIII The 
dashed lines indicate the Icr statistical error bars associated with the chirally extrapolated value 
of the condensate. 
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FIG. 21: Four dimensional norms A^(s) associated with the first eight eigenfunctions of Dh, mea- 
sured on a typical background gauge field configuration generated on an 8'^ x 8 lattice with f3 = 2.3, 
Lg = 24 and ruf = = 0.02. 
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FIG. 22: Magnitude of the matrix elements {A'jj\T s\Ah) , measured on a typical background gauge 
field configuration generated on a 8^ x 8 lattice with /3 = 2.3, = 24 and m/ = m„ = 0.02. 
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FIG. 23: Lowest 10 eigenvalues of Djj^ on a typical background gauge field configuration generated 
on an S'^ X 8 lattice with /3 = 2.3, Lg = 24 and ruf = 0.02. Solid lines represent results from a fit 
to Eq. M 
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FIG. 24: Distribution of 6m values obtained from a fit to Eq. [53]for an 8^ x 8 lattice with (3 
Ls = 24 and m/ = 0.02. 
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FIG. 25: Scatter plot of |A| verses 5m values obtained from the m^, dependence of A^, which have 
been fit to Eq. [53l Eigenvalues where obtained for an 8^ x 8 lattice with /3 = 2.3, Lg = 24 and 
THf = 0.02. The dashed line represents the curve: A = 5m. 
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FIG. 26: Four dimensional norms A^(s) associated with the first eight eigenfunctions of Dh, mea- 
sured on a typical background gauge field configuration generated on an 8'^ x 8 lattice with f3 = 2.3, 
Ls = 2A,mf = 0.02 and niy = -0.25. 
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FIG. 27: Magnitude of the matrix elements {A'jj\T s\Ah) , measured on a typical background gauge 
field configuration generated on a 8^ x 8 lattice with /? = 2.3, Lg = 24, m/ = 0.02, and m„ = —0.25. 
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FIG. 28: Four dimensional eigenvalue density as a function of lambda for 8^ x 8 lattices with 
P = 2.3, two values of Lg and mf = 0.02. Results were obtained from the mj dependence of A'jj, 
which have been fit to Eq. [53l The dashed vertical lines represent the maximum value of A for 
which the distribution is reliably computed. The dashed horizontal line corresponds to the chiral 
limit extrapolation value of the gluino condensate. 
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FIG. 29: Integrated condensate obtained from Eq. [SHwith = ruf for eigenvalues in the range 
— A. Solid curves indicates the integrated condensate as a function of A for 8^ x 8 lattices with 
f3 = 2.3, three values of Lg and mj = 0.02 The corresponding horizontal dashed lines indicate the 
value of the gluino condensate obtained from Eq. |30l as well as the = oo extrapolated value of 
the condensate obtained in Table 1X1 for the fit range 24 — 48. 
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FIG. 30: Distribution of values for an 8^ x 8 lattice with (3 = 2.3, Lg = 24 and mj = 0.02. 
Results were obtained from the niy dependence of A^^, which have been fit to Eq. [53l 
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